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Abstract — We give numerically tractable, explicit integral ex- 
pressions for the distribution of the signal-to-interference-and- 
noise-ratio (SINR) experienced by a typical user in the down- 
link channel from the fc-th strongest bases stations of a cellular 
network modelled by Poisson point process on the plane. Our 
model comprises power-low path-loss model with arbitrarily 
distributed shadowing, with and without Rayleigh fading. Our 
results are valid in the whole domain of SINR, in particular for 
SINR < 1, where one observes multiple coverage. In this latter 
aspect our paper complements previous studies reported in [1]. 

Index Terms — Wireless cellular networks, Poisson process, 
shadowing, fading, SINR, multiple coverage, symmetric sums. 

I. Introduction 

Shannon's theory and its modern extensions quantify the 
quality of communications channels (ergodic capacity, finite 
block errors, error exponents, etc.) in a probabilistic manner 
by considering averages over codewords and channel char- 
acteristics (noise, fading, etc.). For communication networks 
with many channels, it has been recently suggested to use a 
stochastic geometric approach [2] consisting in taking spatial 
averages over node (emitter, receiver) locations. Establishing 
clear connections between stochastic-geometric averages and 
basic information-theoretic notions may be difficult (cf e.g. [3]) 
but this approach for wireless networks has recently attracted 
a lot of attention. In particular, the basic channel characteristic 
discovered in information theory — the signal-to-noise-ratio 
— is now being re-defined and studied in many geometric 
contexts. This paper contributes to this line of thought by 
considering the distribution of the signal-to-interference-and- 
noise-ratio (SINR) of a typical user on down-link channels 
from different bases stations of a cellular network modelled 
by Poisson point process on the plane. In particular, it com- 
plements [1], by providing explicit characterization of low (less 
than one) values of SINR. Current cellular-network technology 
allows for effective use of such SINR regimes, whence our 
motivation comes. 

Cellular network models based on the Poisson point process 
have been shown to give tractable and accurate solutions [4], 
with the Poisson assumption being justified by representing 
highly irregular base station deployments in urban areas [5] 
or mimicking strong log-normal shadowing [6], or both. 
Knowledge of the distribution of SINR allows to calculate 
key performance indicators of cellular networks, e.g. spectral 
efficiency [7] 1 or energy efficiency [6, 9]. It can also be (less 
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immediately) related to user-level quality-of-service metrics. 

In studying the SINR of a given user with respect to the 
best base station (offering him the smallest path-loss), the main 
difficulty for small values of SINR is the necessity of taking 
into account multiple coverage. Indeed, simple algebra shows 
that there can be at most one base station offering a given 
user SINR > 1, and hence the probability of having at least 
one station covering at this level boils down to the sum of 
probabilities of SINR-coverage over all base stations. This is no 
longer the case for SINR < 1, where one indeed needs to study 
probabilities of simultaneous coverage by several base stations. 
We express these probabilities via the so-called symmetric sums 
and relate them to the appropriate partitioning of the SINR 
domain for SINR < 1, which are the main ideas behind this 
paper. 

Our path-loss model consists of the deterministic power- 
law path-loss function and independent, arbitrarily distributed 
shadowing. As previously observed [6, 10], any characteristic 
involving only the sequence of path-loss values experienced 
by a given user from all base stations entirely depend on the 
distribution of the shadowing via its moment of order 2//3, 
where f3 is the path-loss exponent. We also study the impact 
of Rayleigh fading on the SINR coverage by incorporating into 
the path-loss model additional independent random variables 
with exponential distributions. Assuming that fading affects the 
SINR coverage condition but not the choice of the serving base 
station, we observe that it worsens the coverage particularly at 
small values of SINR. 

Related work 

The SINR coverage in a K-tiei network was studied in [1] 
for SINR > 1. Two different approaches to express the 
distribution of the SINR in its whole domain, both involving 
inversion of Laplace transforms, were presented in [6, 11]. Our 
expressions, which involve two key families of integrals (over 
the positive real line and a hyper-cube respectively) are much 
more tractable in numerical evaluation. Our model with fading 
was recently examined in [12] under slightly more general 
assumptions. We revisit it in order to present more closed- 
form expressions under our specific assumptions, and compare 
the coverage probabilities obtained in it to these in the model 
without fading. 2 

2 Which of these two models is more appropriate depends on the context. For 
example if users are motionless, then the fading should be taken into account, 
while highly mobile users "see" channel characteristics averaged over fading, 
cf [13, Proposition 3.1]. 



II. Model description 

On R 2 , we model the base stations with a homogeneous 
Poisson point process $ with density A. Given $, let {S x } xe $ 
be a collection of independent and identically distributed 
random variables, that represent the shadowing experienced 
between station x g $ and a typical user located, without loss 
of generality, at the origin. Let S denote a generic shadowing 
variable. The distribution of S is arbitrary except for a technical 
assumption E[S*?] < oo and a conventional assumption that 
~E[S] = 1 that we make without loss of generality. 

A. SINR multi-coverage 

We define the SINR of the typical user with respect to the 
station x <E $ by 

S x /i{\x\) 



SINR(ir) := 



(1) 



w + i-s x /i{\x\y 

where the constant W is the noise power, / = X^e* S x /((\x\) 
is the total power received from the entire network, and the 
path-loss function is 

l(\x\) = {K\x\f, (2) 

with constants K > and /3 > 2. In this paper we are 
interested in the distribution of the coverage number of the 
typical user defined as the number of base stations that the 
typical user can connect to at the SINR level T, namely 

N{T) = J2 1 [SINR(z) > T] . (3) 

16$ 

The probability of the typical user being covered by at least k 
base stations, which we call k-coverage probability, is 

pW(T) =P{Af(T) > k}. (4) 

In particular, the coverage probability of the typical user is ( 3 ) 
P e (T) :=P C (1) (T). 

Since the function x/(A — x) — Aj (A — x) — l is increasing 
in x it is easy to see that Pc k \T), as a function of T, is 
the tail-distribution function of the SINR experienced by the 
typical user with respect to the base station offering to him the 
kth smallest path-loss Y k : (Y k y 1 /(W + I - (Yfc) -1 ), where 
Y\ < Y-z < . . . is the process of order statistics of {l(\x\)/S x : 
x € $}. In particular, P C (T) is the tail-distribution function of 
the SINR with respect to the smallest path-loss base station. 

Related quantities of interest include also the expected 
converge number E[7V(T)] = J2n=o nP { J ^( T ) = n } and 
its probability-generating function G(z) = E[^^^]. 

B. Adding fading to the model 

In this extension of the previous model we assume that 
the path-loss of each base station x £ $ is further modified 
by a random fading variable F x and equal to l(\x\)/(S x F x ), 
where given <&, {F a } ae $ is a collection of independent and 
identically distributed random variables, independent of shad- 
owing {S x } x( zq,. In this paper we will assume Rayleigh fading, 
i.e., that the generic fading variable F is exponential, with 
E[P X ] = 1. A key assumption is that fading perturbs the 
SINR coverage condition but not the choice of the serving 

3 This notation is similar to that of [1], who use /3 instead of T. 



base station 4 . In consequence, the coverage probability under 
fading (with respect to the smallest fading-averaged-path-loss 
base station) is defined as 

(Y 1 )-'F 



Pc(T) = P 



> T 



(5) 



W + I — (Y 1 )~ 1 F 
where, recall, Yj is the smallest path-loss received by the 
typical user in the model without fading. 

III. Preliminary Observations 

A. Invariance with respect to the shadowing distribution 

Lemma 1: [Cf [6, 10]] The fading-averaged-path-loss pro- 
cess (Y n : n > 1), considered as a point process on the positive 
half-line IR + is a non-homogeneous Poisson point process with 
intensity measure A ([0, t)) = at^ where 

_ AttE[5?] 



K 2 



(6) 



Consequently, the distribution of (Y n : n > 1), and hence the 
functions Pc(T) and P C {T), depend on the model parameters 
(including the shadowing S distribution) only though the noise 
level W, path-loss exponent f3 and the constant a. 5 

B. Symmetric sum representation 

For any given T and n > 1 define the n th symmetric sum 
5„(T):=E[ Y, P{SINR(aJi)>r,i = l,...,n|#}], 

distinct 

(7) 

where P{...|$} denotes the conditional probability given 
$ (it involves random shadowing marks). One defines also 
Sq(T) = 1. We have the following identities related to the 
famous inclusion-exclusion principle (cf e.g [14, IV.5 and IV.3] 
for (8) and (9), respectively). 6 
Lemma 2: We have for k > 1 

oo / _ i \ 

p? ] (T) = E(- 1 ) n " fc (fr_iK( T )' (8) 



P{Af(T) = k} 



EM)' 



-k 



n—k 



S n (T). 



(9) 



n=Q 



(z - l) n S n (T) , z e [0, 1] (10) 



E[Af(T)} = 5i(T). (11) 

Our goal in Section IV-B will be to evaluate symmetric sums 
S n (T), which will allow us to express easily our quantities 
of interest appearing in the right-hand side of the above 

4 In other words, the user compares the received signals averaged over fading 
effects, which is justified by short time and space coherence properties of the 
(multipath) fading. 

5 This means that evaluating our quantities of interest in our model with a 
general distribution of shadowing S and some value of the constant K, one 
can equivalently and for mathematical convenience assume some particular 
distribution S of shadowing, e.g. exponential or constant, and K = 1 provided 
one replaces A by \{K / K) 2 E{SP]/E[SP 1 in the obtained formula. We will 
use these two representations of the general models when evaluating P c (T) 
and P C (T), respectively, in Sections IV-B and IV-C. 

6 A general relation between the distribution of M and the symmetric sums 
is given by the Schuette-Nesbitt formula, often used in insurance mathematics, 
cf [15]. 



expressions. Before doing this, in the following section we 
explain that the (apparently infinite) summations presented 
above boil down to finite sums, as for any given T we have 
Sn(T) = for n large enough. 

C. Partition of the T -domain 

For real x denote by \x~\ the ceiling of x (the smallest integer 
not less than x). 

Lemma 3: For n > 1, S n (T) = whenever T > l/(n— 1). 
In other words, one can replace oo by |"1/T] in the sums in 
expressions given in Lemma 2. 

Proof: This is a consequence of the well known constraint 
of the SINR cell intersection. Indeed, if the SINR of a given 
user with respect to n distinct stations is to be larger than 
T, then n,T/(l + T) < 1 (cf [16, Proposition 6.2], with the 
strict inequality holding whenever there are other (interfering) 
stations or external noise, which is the case in our model. 
Consequently, for T > l/(n — 1) all the terms (probabilities) 
in the mth symmetric sum S m (T) in (7) are null for every 
m> n. ■ 

IV. Main results 

A. Key integrals 

We now introduce two families of functions which will allow 
us to express S n (T) and, in consequence, the multi-coverage 
characteristics in the model without fading. For x > define 

2" f n °° u 2n-l e - M 2 -„^r(l-2/0)-^ 2 dw 
Zn,/»(a) = " ^TTT^n^: Tsl (12) 



j 8 n ^ 1 (C / (/8)) n (n- 1)! 



where 



2tt 



C ^ /3sin(27r//3) 
Remark 4: We have 

Zn,/3(0) = 



= r(l - 2//3)r(l + 2/0). (13) 



The second family of functions are integrals over the hyper- 
cube. For x > define 

^(W)- 1(1 _^ )2//3 

Jn,p{x) = \ — —. d,Vi . . . dv n -! 

7[0 ' 1] ^ 1 Uix + TH) 

i=l 

(15) 

where m := (1 - v. t ) n*=H-i v k- 

Remark 5: For J%, a closed-form solution exists 

J%A X ) (16) 
g(2/g + 1, 2//3 + l) 2 Fi(l, 2/1 + 1; 2(2//? + 1); -l/x) 

a; 

where B is the beta function [17, eq. 5.12.1] and 2-Fi(&, c; z) 
is the hypergeometric function given by [17, eq. 15.11] (whose 
integral representation follows from eq. 15.1.2 and 15.6.1 
therein). 7 

7 We note that the form of J n ,p is similar to integral representations of the 
generalized hypergeometric function and a related integral generalization [18]. 
A closed-form solution of tT n p{x) may exist, but that is left as a future task. 
For low and intermediate n, regular numerical and Monte Carlo methods work 
well and give results in a matter of seconds on a standard PC machine; cf [19]. 
For high n, analysis of the kernel of J n may lead to judiciously choosing 
suitable lattice rules, thus allowing for relatively fast integration [20]. 



B. Results for the model without fading 
For 0<T<l/(n— 1) define 

T 



, , lVr -- d7) 
1 — (n — 1)1 

We now present the key result for the model without fading, 
which gives an explicit expression for the symmetric sums 

S n {T). 

Theorem 6: Assume shadowing moment condition 
E(5 2 /' 3 ) < oo. Then 

S n {T) = T- 2n ^l n ^(Wa-^ 2 )J n AT n ) (18) 

for < T < l/(n — 1) and S n (T) — otherwise, where a is 
given by (6) and T n by (17). 

Theorem 6 in conjunction of Lemma 2 and Lemma 3 give 
us in particular the following expression for the fc-coverage 
probability: 

Corollary 7: Under the assumptions of Theorem 6 



fl/T] 



P^'(T)= ^ (-ir- k (rl)T~ 2 "^X n AWa-^)J n MT n ), 
n—k 

Proof of Theorem 6: Assume < T < l/(n — 1) 
(otherwise the result follows from Lemma 3). Following the 
remark in Footnote 5 we will first evaluate S n (T) assuming 
exponential distribution of S and K = 1, and then bring back 
the general assumptions appropriately rescaling the Poisson 
intensity A. By the (higher order) Campbell's formula and the 
Slivnyak's theorem (see [16, (9.10) and (9.16)]) and a simple 
algebraic manipulation 

Sn(T) = y P{ f] (SINR'fo) > T') }ndn . . . r n dr n 



(19) 



(14) where SINR' (n ) := 



(W+I+ j: U s M "d r := T/(l+T) 
with / is as in (1) and Si exponential (mean-1) variables, 
mutually independent and independent of /. Moreover, the 
event whose probability is calculated in (19) is equivalent to 



> T 



'} 



(20) 



l w + i + YTj^Sj/t^) 

For integer i g denote Ei :— Si/£(ri). By our 

previous assumption Ei are independent exponential vari- 
ables with means 1/pi = l/£(ri), respectively. Let Em '■= 
min(i? 1 , E 2l ■ ■ ■ , E n ). Note that Em is exponential variable 
with mean 1/pm = l/Q~^rLi Mi)- Moreover, define the 
random variable D := 2~27=i &i — tiEm- By the memory- 
less property of the exponential distribution note the random 
variable D is independent of Em and has a mixed exponential 
distribution characterized by its Laplace transform 



l 



(21) 



Using the new random variables we can express the event (20) 
as {E M > T n (W + I + D) } where T n is given by (17). 
Consequently, the probability P{...} calculated in (19) is 
equal to C w (^MT n )Ci(pMT n )C D (pMT n ), which is a prod- 
uct of three Laplace transforms. The first transform is simply 



Cw{0 = e > me second can be shown ([16, equation 
2.25]) to be £/(£) = e -M 2/l, *c'([i)/K 2 while the last one 
is given in (21). After substituting the explicit path-gain 
function (2), noting that there is some symmetry in the in- 
tegration variables r i7 changing the integration variables s ; := 
niXTn^TrC'iP)) 1 / 2 and replacing A by a/(7iT(l + 2/0)) 
to revoke the exponential shadowing assumption and bring 
back the general distribution of shadowing and constant K 
(cf Footnote 5) one obtains 



S n (T) 



T„ 2n//3 (C'(/3))«(n-l)!io Jo 



oo />oo 



(22) 



*=i { Y{ S^e-WiaT {1-2/ P)y 

\i=l 



±4 



nK + r„E»£] 

i=2 fc=l 



dsi . . . ds r , 



(23) 



A substitution of n-dimensional spherical-like variables (de- 
tailed in the appendix, Section A) completes the proof. ■ 

C. Effects of Rayleigh fading 

We now consider the model with fading; cf Section II-B. 
Theorem 8: The coverage probability under fading ( defined 
in (5) is equal to 



Pc(T) 



2 

P 

x exp 



-tTWa^ 13 ' 2 -t$ 



2 Tt 2 /P 2 F 1 (1,1 

"p 



2/0; 2- 2/0; -T) 



(1-2/0) 



(24) 



dt. 



where, again, 2F\(a,b; c; z) is the hypergeometric function 
mentioned in Remark 5 above. 

Remark 9: The expression (24) can be easily evaluated 
numerically. Setting W = yields an analytic solution 



Pc(T) = 



1 



2 2^1 (1,1 





2/0; 2- 2/0; -T) 



(1-2/0) 



T 



(25) 



Proof of Theorem 8: We user the path-loss process 
representation {Y n } defined in Section II-A (which does ac- 
count for arbitrary general shadowing, but not for fading), 
which we enrich by independent exponential marking F n 
representing Rayleigh fading. By Lemma 1 {(Y n ,F n ) : n > 1} 
is independently marked Poisson point process of intensity 
A( ). Using this representation we can express the coverage 
probability of (5) as follows 



P C {T) = j P{ F 1 > sT(W + / (S)00) ) }f Yl (s) 



ds 



£ w (sT)£( S!OC )(sT)f Yl (s) ds , 



(26) 



where /^(s) is the probability density of Yi, known to be 
(due to Poissonianity of {Y n }) fy^ds) = -^e~ A(s) = 
2a/(3s 2 / l3 ~ 1 and where I( s ,<x) i s trie random variable repre- 
senting conditional interference (accounting for shadowing and 
fading) given Y\ — s. Again, it is well known that I/ S ,oo) i s 
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Fig, 1: Distribution function of SINR from the strongest base station with and without 
fading validated by the Laplace inversion method and simulation. 

equal in distribution to J^y >s Fn/Y n and has the Laplace 
transform which can be explicitly evaluated as follows 



exp 



[l-£ F ({/v)]A{dv)\ (27) 
2a £s 2 //3 2 F 1 (l,l-2/fr2-2/ft-€/s 



(1-2/0) 



(28) 

Plugging into (26) and substituting t = so! 5 ! 2 completes the 
proof. ■ 

V. Numerical illustrations 

We use MATLAB implementation [19] for all our calcula- 
tions. We set /3 = 3.8 and K = 6910 km -1 (which corresponds 
to the COST Walfisch-Ikegami model for urban environment). 
The shadowing is modeled by a log-normal random variable of 
expectation 1 and logarithmic standard deviation lOdB (cf [6]) 
which makes F,(S 2 ^) = 0.516. We assume noise power 
— 96dBm normalized by the base station power 62.2dBm which 
makes W = 10 - 15,82 . We consider two values for the density 



of base stations: A = 4.619km , which corresponds to 
a "urban" network deployment and A 



0.144km" 2 for a 



"suburban" one. Figure 1 shows the distribution function of 
SINR from the strongest base station for both scenarios. We 
validate our approach by showing that the obtained results 
coincide with those of simulation and a Laplace inversion 
method developed in [6], with the latter approach being less 
numerically stable and much more time-consuming. We also 
plot the distribution of SIR in both scenarios (i.e. assuming 
W = 0). Both SIR curves coincide with that of the SINR in 
urban area, thus showing that for urban density of stations the 
network is interference-limited, while for suburban density the 
impact of noise is non-negligible. Finally, we provide curves 
regarding the model with fading (i.e. 1 — P C (T)). We observe 
that the impact of fading is non-negligible in both the urban 
and suburban scenario, and stochastically decreases SINR (the 
respective distribution functions are larger). 

VI. Conclusion 

Cellular network models based on the Poisson point process 
allow for analytic expression for many important characteris- 



tics. Complementing previous studies, in this paper we give 
tractable, integral expressions (no Laplace transform inversion 
needed) for the distribution of the SINR experienced by a typ- 
ical user in the down-link channel from the k-th strongest base 
station. Our model comprises power-low path-loss model with 
arbitrarily distributed shadowing, with and without Rayleigh 
fading. 

Appendix 

A. Remaining proof of Theorem 6 

We introduce a change of variables inspired by the n- 
dimensional spherical coordinates (for example, see [21, eq. 
(1.3)]) 



si 

S2 
«3 



u[sm V\ sin t>2 ■ ■ ■ sin tV„_i 



. sin n -\ 



[2//5 
12//3 



|2//3 



it [cos v i sm 02 ■ ■ 
u[cos 02 sin 03 . . . sin n -\ 

s n := u[cos0 n - 1 } 2/ ' 3 . 

Observe that J2ti 4 = and U7=i «i = «" IlLi %f //? , 
where g f = . . . , 0„_i) := (s i /u)' 3/2 . When /3 = 2 our 

system of coordinates boils down to the regular n-dimensional 
spherical coordinates, whose Jacobian is J(u, 0\, . . . , 6 n -i) = 
u"- 1 n^i 1 sin^" 1 6i\ cf [21, eq. (1.5)]). By induction (or 
determinant properties and the chain rule) our coordinate 
system has the corresponding Jacobian 



J(u,6i, . 

2 

1 



, 0n-l) 
-1 

J(U,01, . 



(29) 



■ ,0 n ) 



'n-l 

n 



-i 2//3-1 



sin 



.; COS < 



(30) 

which is clearly postive over the integration domain of interest. 
Denote z := W(aT(l - 2/P))~^/ 2 . The integral in (23) 
becomes 



n-l 

\~[ sin 2 0i cos 0i 

i=l 



[0,7r/2]»- 
X J(u, 0\, 



i=2 



2(/3+l)//3 



n ) du d0\ 

oo 

2n-l. 



(I 



n [s 



, COS ( 



1 4//3+1 



sine 



[0,tt/2]* 



nk 2 



-d6»i . . . d(9„ 



The substitution = sin 6i makes the second integral (over 
the hypercube) equal to 2 1 ~ n j7„ j/ g(T„), which, after defining 
rji, completes the proof in view of (23) and (17). 
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